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ABSTRACT 



r-^ , Aims. We discuss the current X-ray view of the Crab Nebula and Pulsar, summarising our analysis of observations of the source with the 

^^ EPIC-pn camera on board the XMM-Newton observatory. Different modes of EPIC-pn were combined in order to yield a complete scenario 

C^ of the spectral properties of the Crab resolved in space and time (pulse phase). In addition we give a description of the special EPIC-pn Burst 

^^^ mode and guidance for data reduction in that mode. 

^^ Methods. We analysed spectra for the nebula and pulsar separately in the 0.6-12.0 keV energy band. All data were processed with the SAS 

^5 ' 6.0.0 XMM-Newton Scientific Analysis System package; models were fitted to the data with XSPEC 1 1 . The high time resolution of EPIC-pn 

f^ ^ in its Burst mode (7yus) was used for a phase resolved analysis of the pulsar spectrum, after determination of the period with epoch folding 

I techniques. Data from the Small Window mode were processed and corrected for pile-up allowing for spectroscopy simultaneously resolved 

O in space and time. 

-4-^ • Results. The spatial variation of the spectrum over the entire region of the Crab shows a gradual spectral softening from the inner pulsar region 

^ ■ to the outer nebula region with a variation in photon index, F, from 2.0 to 2.4. Pulse phase resolved spectroscopy of the Crab Pulsar reveals a 

. . , phase dependent modulation of the photon index in form of a significant hardening of the spectrum in the inter-peak phase from F = 1 .7 during 

^ ■ the pulse peak to F = 1.5. 

X' 

d Key words, stars: neutron stars - pulsars: individual: PSR 0531-1-21 - supernova remnants: Crab - X-rays:stars, - instruments: EPIC-pn, - 
data analysis: Burst mode 

1. Introduction lished from RXTE ("Pravdo. Angelini. & H ardind. fl997l) and 

from Beppo-SAX ( Massaro et al.^ .200 0). Both measurements 

Since the discovery of the Crab Pulsar in 1968 showed a significant hardening of the spectrum in the inter- 

j Staelin&Reifensteint |I968V the Crab has been one of p^^^ ^^^.^^ ^^^ ^^^ Beppo-SAX both were unable 

the best studied objects in the sky and it remains one of the ^^ ^^^^j^^ ^j^^ ^^^^ ^p^^^^jl^ ^^^^^^ ^^^^^^^ observations 

brightest X-ray sources regiilarly observed. With an X-ray have provid ed spatially resolved spectroscopy of the Crab 

luminosity of ~5 lO^^ergs"', the Crab emits a significant JMori et al .-.'2004). These observations, however, suffered from 

fraction of its energy output in the X-ray band (the total pUe-up which had to be especially corrected for Similar correc- 

luminosity of the Crab is ~2 10 ergs ). As a standard candle .. ., ■ ,. .f ., f ., cnrr^ c n 

^ ° ' tion methods are discussed for the case of the EPIC-pn Small 

for instrument cahbration, it has been repeatedly studied window mode in Sect. E3| 
by many astronomy missions. The 33 ms Crab Pulsar has 
been observed in almost every band of the electromagnetic 
spectrum. Its pulse profile exhibits a double peaked structure With t he European photon Imaging Camera (EPIC; 



with a phase separation of 0.4 between the first and the second IStriideret a l.. 2001; Tur ner et all l200ll) in its Burst mode, 

pulse. XMM-Newton provides the unique possibility to measure the 

Pulse phase resolved X-ray spectroscopy of the Crab was Crab spectrum without pile-up when combining one dimen- 

perfo rmed for the first time with OSO 8 ( Pravdo & Serlemitsos, sional spatial resolution with a high time resolution (7 //s), al- 

Il98lh . Analyses with a finer phase resolution were pub- lowing for spectroscopy simultaneously resolved in space and 

time (phase). In addition, we derived two-dimensionally re- 

Send ojfprint requests to: M.G.F. Kirsch, mkirsch@sciops.esa.int solved spectra using bootstrapping techniques and Monte Carlo 
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simulations to correct XMM-Newton EPIC-pn Small Window 
mode data. 

This paper is organised as follows. In Sect. |2l we give a 
description of our observations and the XMM-Newton instru- 
ments, followed by some technical comments on data analysis 
of Burst mode observations in Sect. |3l Our results from timing 
and spectral analysis of the Crab are presented in Sect.|4l 

2. Observations performed by XMM-Newton 

XMM-Newton JJansen et all I2OOII) was launched in 1999 
December with an Ariane 5 rocket from French Guyana. It op- 
erates six instruments in parallel on its 48 hour highly ellipti- 
cal orbit: three Wolter type 1 telescopes, with 58 nested mir- 
ror shells each, focus X-ray photons onto the five X-ray instru- 
ments of the EPIC ( Struder et al. , 200 1 ; Turner et al. , 200 1 ) and 
the R eflecting Grating Spectrometers (RGS: Iden Herder et all 
I2OOII) . In addition, a 30 cm Ritchey Chretien optical tele- 
scope, flie_0£tical Monitor, is used for optical observations 
(OM; iMasonetall 1200 ih . EPIC consists of three cameras: 
The two EPIC-MOS cameras use Metal-Oxide Semiconductor 
CCDs as X-ray detectors, while the EPIC-pn camera is 
equipped with a pn- CCD. Both have been especially developed 
for X MM-Newton ('Pf effermann et all [l999t iMeidinger et all 
[1999 : Turner et al., 2001^ 

EPIC provides spatially resolved spectroscopy over a field- 
of-view of 30' with moderate energy resolution. The EPIC 
cameras can be operated in different observational modes re- 
lated to different readout procedures. Detailed descriptions of 
the various readout modes of EPIC-pn and their limitations 
are given by lKendziorra et al. ( 1999), Kuster et al. ( 1999) and 
lEhle et all ( 120031) . The EPIC-pn camera, which provides the 
highest time resolution (Timing mode: 30fj.s, Burst mode: 
7jus) and moderate energy resolution (E/dE = 10-50) in 
the 0.2-15 keV energy band is the ideal instrument to ob- 
serve the Crab. In Burst mode the pile-up limit for a point 
source is 60000 counts s"', corresponding to a maximum flux 
of 6.3 Crab. The special readout of the Burst mode, how- 
ever, leads to a loss of spatial resolution in the shift direction. 
Moreover, the lifetime in Burst mode is only 3%, limiting the 
mode to observations of very bright sources such as the Crab 
or very bright transients. 

The Crab was observed with the EPIC-pn camera in 
Burst mode during XMM-Newton's revolutions 234 and 411. 
Observations were also performed in the EPIC-pn Timing and 
Small Window modes, however the high count rate of the 
source implies severe pile-up effects which distort the spectrum 
and require special treatment in the analysis. The Timing mode 
data have not been used here, since they do not contain addi- 
tional information compared to the Burst mode data. The Small 
Window mode data, however, are used despite of the pile-up 
distortion, as they provide a two-dimensional spatial resolution 
of the source. 

Table ^ summarises the observations of the Crab with 
EPIC-pn which were included in our analysis, listing opera- 
tional modes and total exposure times for each individual ob- 
servation. The exposure times were corrected for data losses 
due to telemetry constraints and detector deadtime. 



Table 1. Log of observations of the Crab used in this paper. 



OBSID 


Rev. 


Exp. 

[s] 


Mode 


Filter 


Position 

Angle[deg]' 


0122330701 


56 


166.5 


SW 


T 


269.6 


0135730701 


234 


298.1 


BU 


T 


269.3 


0153750201 


411 


105.1 


BU 


M 


267.3 


0153750301 


411 


181.7 


BU 


M 


267.3 


0153750501 


411 


201.6 


BU 


M 


267.3 



BU: Burst, SW: Small Window, T: Thick, M: Medium, Exp.: Exposure 
time, taking also the life time of the different modes into account. 
* Rotation between the spacecraft A'-Z-plane and the plane defined by 
the spacecraft X-axis and celestial North. 
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Fig. 1. The special readout in Burst mode, which gives the pos- 
sibility to observe very bright sources with negligible pile-up. 

3. Data analysis for EPIC-pn Burst mode 

In the following section we describe the special techniques of 
data reduction needed for the EPIC-pn Burst mode in order to 
be able to derive spectra and pulse profiles of the Crab. 

The data sets were processed using the XMM-Newton 
Scientific Analysis System, SAS 6.0.0. Event times were cor- 
rected to the solar system barycentre using the SAS tool 
barycen. Single and double events were considered for the 
analysis while FLAG = was chosen to exclude border pixel 
events, for which the pattern type is uncertain and thus, the total 
energy is only known with lower precision. 

The Burst mode operates a special readout similar to a tape 
recorder. Within 14.4yus, 200 lines are fast-shifted while accu- 
mulating information from the source. Similar to the Timing 
mode, this leads to the loss of spatial resolution in the shift di- 
rection. The stored information is then read out in a standard 
fashion, ignoring the last 20 lines of the CCD because of con- 
tamination by source photons arriving during the readout. The 
CCD is then erased by a fast shift of 200 lines and immedi- 
ately afterwards the next Burst readout cycle starts. Fig. [2 il- 
lustrates this readout procedure. In the following we use SAS 
terminology for the images, i.e., RAWY represents the CCD 
line numbers, i.e., the position on the CCD in shift direction 
from 1-200 and RAWX represents the CCD column numbers, 
i.e., the position perpendicular to the shift direction from 1-64. 

Note that events with RAWY > 180 cannot be used for 
any kind of data analysis, since they contain severely piled- 
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Fig. 2. Variation of spectral parameters with extraction region 
in EPIC-pn Burst mode. Spectra have been extracted in a sUd- 
ing box of 13 RAWY pixels in steps of 10 RAWY pixels. An 
absorbed power law fixing A^h = 2.76 10^' cm"^ (see Table|5J 
derived from the nebula region of the Crab was fitted to the 
data. In addition we show the variation of the spectral parame- 
ters extracting spectra from a growing box in order to visualise 
the effect as averaged over columns. Spectra should only be 
extracted for RAWY < 160 (for determination of the photon 
index) and RAWY < 140 (for determination of the normalisa- 
tion). The background was also taken from a sliding box that 
contains in first approximation no source photons. Note that the 
dashed red line without symbols is just to guide the eye with re- 
spect to the normalisation parameter variation of the growing 
box (dashed red one with symbols). 



up data as accumulated from the point-spread function (PSF) 
wings during the slow readout phase of a Burst cycle. (The 
HEW of XMM-Newton's EPlC-pn PSF is 15" and one EPIC- 
pn pixel has a size of 4.128x4.128n"). These data are there- 
fore not transmitted to ground at all. We quantified the effect 
of pile-up in CCD lines 1 < RAWY < 180 for the measured 
spectrum of the Crab by extracting spectra with a sliding box of 
13 RAWY pixels length, shifted in steps of 10 RAWY pixels. 
An absorbed power law fixing A^h at 2.76 10^' cm"^, as derived 
from the nebula region of the Crab (Table |2ji, was fitted to the 
data. The fit results show that the power law index is only af- 
fected for RAWY > 160 while the normalisation seems to be 
already significantly affected as of RAWY ~ 140. Therefore 
we recommend for general data analysis in the EPIC-pn Burst 
mode to extract spectra only from RAWY-regions < 160 for 
determination of the spectral slope and RAWY < 140 for de- 
termination of the normalisation. 

The effects due to the special readout in Burst mode on 
both the spectral shape and the normalisation of the flux are 
shown in Fig |2l For our spectral analysis we excluded events 
with RAWY > 142. 

The Burst mode in EPIC-pn provides spatial resolution 
only perpendicular to the CCD readout direction, i.e., in the 
RAWX direction, which depends on the spacecraft's position 
angle with respect to the source. To visualise the orientation of 
the Crab Nebula with respect to the readout direction, we com- 




Fig. 3. Upper panel: Intensity image of the Crab in EPIC-pn 
Small Window mode (rev. 56). Lower panel: Crab image in 
Burst mode. All events in the full energy range (0.2-15 keV) 
were used. The solid box indicates the pulsar extraction re- 
gion, while the dashed and dotted boxes are used for nebula 
and background data extraction respectively. 



pare in Fig.|3la Burst mode image of the Crab (revolution 234) 
to an image obtained in Small Window mode (revolution 56) 
which was observed under nearly the same position angle (see 
Table nj. The additional observations in revolution 411 were 
done with the slight difference of 2? 3 in the position angle. 
This variation, however, is still negligible for the purpose of our 
analysis. Thus, we can directly compare the two-dimensional 
image of the Crab taken in Small Window mode with the one- 
dimensional cut acquired during Burst mode observations. The 
extraction regions to derive spectra of the pulsar, nebula and 
background are also illustrated in Fig.|3] 
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Fig. 4. 0.2-15 keV pulse profile of the Crab 
from the three observations in rev. 411, using 
The profile is shown twice for clarity. 



1.5 

Pulsar obtained 
250 phase bins. 



4. Data analysis results 

4.1. Pulse profile and period 

We determined the period of the Crab pulsar using epoch fold- 
ing techniques. The period P was computed separately for each 
individual observation in revolutions 234 and 411, using ex- 
trapolated values of P supplied by the Jodrell Bank Crab Pulsar 
Monthly Ephemeris (http : //www . jb . man . ac . uk). The best 
pulsar period derived from the data from revolutions 234 and 
411 were 33.52130944(2)ms and 33.5341004590(5) ms for 
the epochs 51988.6405766415 (MJD) and 52340.6825136183 
(MJD) respectively. Thus, the relative deviation of the observed 
pulse period with respect to the most accurate radio data value 
available was AP/P < 10"^. Data were combined for detailed 
analysis from the three observations in revolutions 411. The 
folded 0.2-15 keV pulse profile is shown in Fig.|4l 

4.2. Localisation of the pulsar 

As outlined above, spatial resolution in the RAWY-direction is 
lost when EPIC-pn is operated in Burst mode. To locate the 
pulsar without use of the attitude in such an image, events 
in a running window of 5 columns in RAWX-direction were 
folded with the best-fit period. The pulsed flux for each of 
those windows was calculated subtracting the ofF-pulse region 
(phase 0.568-0.768) as background contribution. The pulsed 
flux reached its maximum, indicating the location of the pul- 
sar, when the window was centred on column 37. In a next 
step, the pulsar spectrum was extracted taking a larger window 
of seven columns width (28'.'896) centred at column 37. The 
spatial variation of the total flux and the pulsed flux is shown 
in Fig.|5] The peaks of the distributions of the total flux and the 
pulsed flux are not co-aligned, indicating that the pulsar is not 
located in the centre of the nebula. The dashed vertical lines 
show the region chosen for the extraction of the pulsar spectra. 
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Fig. 5. Variation of the 0.5-15 keV total flux (continuous lines) 
and the pulsed flux (dashed-dotted lines) over the entire Crab. 
Note that the flux units are per 7 columns window and renor- 
malised to the same extraction area, using a normalisation fac- 
tor of 180/140. For the pulsed flux we subtracted the renor- 
malised mean flux from the off'-pulse region in phase 0.568- 
0.768. The two curves shown are for Rev. 411 and Rev. 234, 
respectively. The difference in flux is caused by the use of dif- 
ferent filters. 




channel energy (keV) 

Fig. 6. Spectra of the overall Crab region covering the 0.6- 
9.0 keV energy band. Data points marked in black and red 
correspond to revolutions 411 and 234 respectively. The back- 
ground contributions were subtracted from the data as listed in 
Table|2l 



4.3. Spectra of the Crab from EPIC-pn in Burst mode 

Spectra from different regions of interest of the Crab were ex- 
tracted from the single Burst mode observation in revolution 
234 and the three merged observations in revolution 411, as 
listed in Table |2 Note that different ways of background sub- 
traction were applied (BGSM). The spectra were fitted with an 
absorbed single power law model spectrum in XSPEC. The key 
spectral fit parameters are also listed in Table|3 Fits and model 
to data ratios are shown in Figs.|6j|9] 
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Fig. 7. Spectra of the pulsar region covering the 0.65-8.5 keV Fig.9. Spectra of the off-pulse region covering 0.6 - 9.0 keV 
energy band. Data points marked in black and red correspond energy band. The data points marked in black and red cor- 
to revolutions 41 1 and 234, respectively. The background con- respond to revolutions 411 and 234, respectively. The back- 
tributions were subtracted from the data as listed in Tabled ground contributions were subtracted from the data as listed in 

TableE 

Table 2. Extracted spectra from different regions of interest and 
spectral parameters. 




Spectrum BGSM* 



A'h 
10^1 cm-2 



xL(dof) 



total Crab 


C-B 


2.60(1) 


2.046(3) 


1.99(3075) 


pulsar 


P-0 


2.9(3) 


1.72(5) 


0.91(75) 


nebula 


N-B 


2.76(4) 


2.263(9) 


1.24(1394) 


off-pulse 


0-B 


2.62(4) 


2.030(9) 


1.23(1296) 



channel energy (keV) 

Fig. 8. Spectra of the nebula region covering the 0.6-9.0 keV 
energy band. Data points marked in black and red correspond 
to revolutions 411 and 234, respectively. The background con- 
tributions were subtracted from the data as listed in Table |2] 

4.3.1 . Spatial variation of tine spectra 

To study the spatial variations of the source in greater depth, 
spectra were extracted from a running window in RAWX- 
direction of five columns width. An absorbed power law model 
was fitted to each spectrum, fixing the equivalent hydrogen col- 
umn density at A^h = 2.76 10^' cm"^ (see Table |2j as derived 
from the nebula region of the Crab. We were able to distinguish 
four typical regions as (a) background (columns 1-13 and 55- 
64) (b) outer nebula (columns 14-22 and 47-54) (c) inner neb- 
ula (columns 23-33 and 41^6), and (d) pulsar (columns 34- 
40). The variation of the spectral index along those extraction 
regions is plotted in Fig. 1101 

The spectrum is hardest in the region of maximum flux and 
pulsed fraction (see Figs.l5landllO>. The total variation of the 
photon index along the Crab is a composed function of the con- 
tributions from pulsar and nebula. While all data points outside 
the pulsar region lie on a parabolic curve, the contribution of 
the pulsar becomes manifest as an additional dip of the photon 



*BGSM: Background Subtraction Method: This column explains 
which region was used to generate the source and which for the back- 
ground spectrum (A-B means: A source; B background spectrum.) 
C: all Crab (columns 22-51), P: pulsar region (columns 34-^0), N: 
nebula region (columns 22-26), B: background region (columns 3-7), 
O: off-pulse region (columns 34-40). 
Note that errors are given in terms of 90 % confidence level. 



index in columns 34-40. Note that Fig.[TO|is a result of spec- 
tral analysis without background subtraction in order to show 
the hardening of the background outside the nebula region. 

4.3.2. Pulse phase spectroscopy 

Events from the pulsar region (RAWY < 142, columns 34- 
40) were sorted into 7 pulse phase intervals. For each phase 
interval we extracted spectra and computed the photon index 
with A^H fixed at the value derived from the off-pulse region, 
A^H = 2.62 10^' cm"^. The reduced ;i'^ values are generally ac- 
ceptable (;if^^^ - 1.2-1.3).Forbackgroundsubtraction, we used 
the off-pulse region (columns 34-^0 with phase 0.568-0.768). 
Plotting the photon index, F, against the phase intervals, we see 
a modulation of F with the pulse phase (Fig.llH. The secondary 
pulse shows a harder spectrum than the main pulse. In the inter- 
mediate pulse region the spectrum is found to be harder than in 
the peak of the pulses. The spectrum of the intermediate pulse 
region is harder than the spectrum of the main pulse by 0.3. 
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^^500 




30 
column 

Fig. 10. Spatial variation of the photon index (stars) over the 
Crab nebula and pulsed flux (solid blue lines) in the energy 
range 0.6-9.0 keV. Note that the flux units are renormalised to 
the same extraction area by a factor 180/140. For the pulsed 
flux we subtracted the renormalised mean flux from the off- 
pulse region (phase 0.568-0.768). Note that since the spectral 
analysis has been performed without background subtraction 
the hardening of the spectral index outside the nebula region 
is influenced by the spectrum of the X-ray background in the 
Crab direction. 
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Fig. 11. Spectral variation with pulse phase of the Crab Pulsar 
Photon index (red crosses) in the energy range 0.6-6.5 keV, 
pulsed flux: blue lines for the two observations in rev 234 and 
411. For this analysis only data of the Burst mode observa- 
tions have been used. The flux has been renormalised to the 
same extraction region by multiplying with a correction factor 
of 180/140. For the pulsed flux we subtracted the renormalised 
mean flux from the off'-pulse region with phase 0.568-0.768 
and took into account that we measure only 77% of the en- 
circeled energy in the extraction region of seven columns. 



The spectrum of the second pulse is harder than the spectrum 
of the main pulse by 0. 1 . 

The analysis of iMassaro et al .' feOOCr) did show, that a two 
component emission model can account for the variation of the 
pulse profile with energy and the behaviour of the spectral in- 



Fig. 12. Two dimensional spectroscopy of the Small Window 
mode data without correction for pile-up. One coloured box 
represents a region of 9 EPIC-pn pixels (153.363 n"), the 
colour represents the spectral index (see colour code scale). 
The count rate contours in Small Window mode are overplot- 
ted. 



dex with phase. However, no complete model for the X-ray 
emission from the Crab Pulsar was provided. 

4.4. Spatially resolved spectroscopy of the Small 
Window mode data 

In Small Window mode, one EPIC-pn CCD is operated in a 
window of 64x64 pixels (Fig.O. Events are read out after an 
integration time of 3.9812 ms. Due to the high source intensity, 
observational data of the Crab taken in Small Window mode is 
expected to be severely affected by pile-up (the limiting count 
rate in Small Window mode before the ons et of pile-up is o f 
the order of 130 counts s ' for a point source. lEhle et al.l2003h . 
However, as outlined above, the Small Window mode provides 
additional spatial information compared to the Burst mode. 
In analogy to the analysis of Burst mode data (Sect. |3l 14. 3> . 
we extracted Small Window mode spectra of the Crab (single 
and double events, FLAG=0) and fitted an absorbed power law 
XSPEC model to the data. The results of our analysis are vi- 
sualised in a colour coded diagram in Fig. [21 The variation of 
the photon index over the Crab is indicated by the colours of 
the boxes, representing extraction regions of 153.363 n" each. 
The observed count rates are overplotted as contours. Note that 
these results are of qualitative nature only and yet have to be 
corrected for pile-up effects. 

In order to estimate the level of data contamination due to 
pile-up we took two different approaches: 

1 . By comparing predicted to measured CCD pattern fractions 
using the SAS task epatplot (Fig. lOt . The change in pattern 
ratios is a sensible indicator of the onset of pattern pile-up 
(Two single photons hitting the detector in the same read- 
out frame in adjacent pixels are misinterpreted as a single 
event of double pattern). Since the lower threshold of the 
EPIC-pn camera was 80 ADU (400 eV) higher than the de- 
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Fig. 13. Pile-up assessment of the Crab observations in Small 
Window mode. We compared the predicted with the measured 
pattern fractions using the &45-task epatplot. epatplot was op- 
erated on a grid of event regions with 7 x7 pixels (834.979 n"). 
The upper panel shows the deficit in single events, the lower the 
increase of double events. Both is an indication for pile-up. 

fault value for the Small Window mode observation we had 
to recalibrate epatplot using pattern fractions of non piled- 
up sources in Small Window mode. This special calibration 
will be provided to the public in the near future via the stan- 
dai-d SAS. 
2. Analysing the loss of events due to pile-up by determining 
the count rate of single columns in Burst mode and compar- 
ing them with the corresponding columns in Small Window 
mode (Fig.O. Note that for this analysis only events with 
RAWY < 1 1 have been used for Burst mode in order to 
avoid any effect of pile-up during readout hampering the 
normalisation (Fig.|2j- Fig.ll4lshows the ratio of the count 
rates in Burst and Small Window mode, where numbers di- 
verging from 1.0 indicate pile-up. 

Both methods confirm that the complete source area of the 
Crab observation in Small Window mode may only be used 
for analysis when correcting the data for pile-up. Pile-up ef- 
fects change both flux and spectral shape of the source data, as 
with increasing count rate the probability increases that several 
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Fig. 14. Count rates and ratio of the count rates in Burst and 
Small Window mode, where numbers diverging from one in- 
dicate pile-up (blue). Black solid: Burst mode count rate, red 
solid: Small Window mode count rate, black dashed: Burst 
mode count rate corrected (a) for pointing difference Burst and 
Small Window observation of 0.5 pixels between, (b) the ex- 
traction region has been reduced to RAWY values from 2-11, 
(c) the lower threshold settings of Small Window mode have 
been also applied to the Burst mode data. Blue solid: uncor- 
rected ratio of Burst and Small Window mode count rate, blue 
dashed: corrected ratio of Burst and Small Window mode count 
rate. 



low-energy events are wrongly detected as single high-energy 
photons. 

In order to quantify the effects of pile-up on the spectral 
shape. Burst mode data and the Small Window mode data have 
been compared in one spatial dimension, considering spectra 
extracted from single CCD columns. This approach simulates 
for the Small Window mode the loss of spatial resolution along 
the shift direction in Burst mode. A direct comparison can thus 
be made for non piled-up data from Burst mode and piled-up 
data from Small Window mode. Fig. [21 shows at the left the 
resulting average photon indices as a function of count rate 
(column) for both Small Window and Burst mode data. The 
piled-up Small Window mode spectra are flatter than the non 
piled-up Burst mode spectra. This comparison gives a first esti- 
mate of the order of pile-up corrections required in general for 
Small Window mode spectra of Crab-like sources, and will be 
used in the following to test our theoretical pile-up corrections 
(Fig. El right). 

This observational estimate, however, can not be used for 
a straightforward correction of the two-dimensional spectral 
analysis, where the observed count rate assumes values up to 
1 count /( s n") in the central region of the Crab. In a more gen- 
eral approach, we performed Monte Carlo simulations, imple- 
menting the pile-up behaviour of the EPIC-pn CCDs. Using 
non-contaminated Burst mode spectra to fix the input spectral 
shape characterised by various photon indices, and varying the 
input flux, we simulated the effect of the Small Window mode 
on the source spectra. Comparing the results from spectral fits 
of Burst input and simulated Small Window output spectra, we 
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Fig. 15. Left: Photon indices as a function of count rate for both Small Window and Burst mode. The piled up Small Window 
mode spectra show lower photon indices than the non piled-up Burst mode spectra for the same count rates. Right: Same data set 
with pile-up corrections for Small Window data. The Small Window count rates are replaced by the correct values from Burst 
mode, the photon indices are corrected according to the function f(r) given by ([Q. Note that the general trend of the data points 
(lower photon index with higher count rate) reflects the spectral behaviour of the Crab as already shown in Fig. 1101 
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Fig. 16. Result of Monte Carlo simulations to estimate the ef- 
fect of pile-up to the measured spectra depending on photon 
index and count rate. 

derived a correction function / for the photon index with re- 
spect to the observed count rate. With the onset of pile-up in 
Small Window mode, the photon index [p.i.^^(r)] decreases 
linearly with the observed count rate r. 
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The negative constant shift (const. ~ 0.03) of the observed 
with respect to the true photon index is caused by the unusually 
high lower threshold of this Small Window mode observation, 
which leads to changes in the fraction of double and single pat- 
terns at low energies and thus to changes of the spectral shape 
in the low energy regime. 

Fig.[^shows the result of our simulations and indicates the 
shape of the correction function for different initial settings. 

The corrections for the photon index have been applied to 
the Small Window mode data. The derived correction function 
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Fig. 17. Distribution of patterns for a merged EPIC-pn Small 
Window mode observation of 620 ks public data. Internal back- 
ground has been subtracted using 125 ks merged Closed-Filter 
observations. Red: singles, blue: doubles, green: triples, light 
blue: quadruples, black: invalid patterns. 

Q is in very good accordance to our independent observational 
quantification of pile-up effects for the one-dimensionally aver- 
aged Small Window data (Fig.^] right). Correction of the two 
dimensional spectroscopy of the Crab yields a final photon in- 
dex of r = 1 .96 (without correction: F - 1.6) for the innermost 
extraction region around the Crab Pulsar 

Note that for that analysis triple and quadruple events can 
be neglected in the case of the EPIC-pn camera since they only 
account for less than 4% of all events, see also Fig. 1171 

5. Conclusion 

Our work provides phase resolved spectroscopy of the mod- 
ulated flux of the Crab Pulsar in combination with spa- 
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tially dependent spectral results of the pulsar and the neb- Willingale, R., et al. 2001, A&A, 365,L212 
ula. We measure a temporal hardening of the spectral index 
in the inter peak region by 0.4. This confirms the results 
achieved by Beppo-SAX (M assaro et al., 2000) and RXTE 
(iPravdo. Angelini. & Hardina 1 19971) measurements. In addi- 
tion, we give some insight into the spatial behaviour of the 
nebula both from one-dimensional EPIC-pn Burst mode and 
two-dimensional EPIC-pn Small Window mode data. We show 
that the pulsar region of the Crab Nebula has by far the hardest 
spectrum getting softer at larger distance from the pulsar. These 
results are in line with the analysi s of XMM-Ne wton EPIC- 
MOS data JWillingale et all l200ll) and Chandra llMori et all 
l2004 . 

For the first time we can verify a simulated pile-up correc- 
tion for XMM-Newton imaging mode data with real data from 
the one dimensionally resolved Burst mode. 

We like to stress that this paper describes in detail how 
XMM-Newton EPIC-pn data from its special Burst mode need 
to be analysed addressing both the excellent prospectives and 
the strength of that mode, but also the caveats that are crucial 
for an adequate treatment of the data and their exploitation. 
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